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Abstract 

In this paper, we study the discrete-time quantum random walks on a line 
subject to decoherence. The convergence of the rescaled position probability 
distribution p(x,t) depends mainly on the spectrum of the superoperator Ckk- 
We show that if 1 is an eigenvalue of the superoperator with multiplicity one 
and there is no other eigenvalue whose modulus equals to 1, then P(-^=, i) con- 
verges to a convex combination of normal distributions. In terms of position 
space, the rescaled probability mass function p t (x,t) = p(y/tx,t), x G Z/y/i, 
converges in distribution to a continuous convex combination of normal distri- 
butions. We give an necessary and sufficient condition for a U(2) decoherent 
quantum walk that satisfies the eigenvalue conditions. We also give a complete 
description of the behavior of quantum walks whose eigenvalues do not sat- 
isfy these assumptions. Specific examples such as the Hadamard walk, walks 
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under real and complex rotations are illustrated. For the 0(2) quantum ran- 
dom walks, an explicit formula is provided for the scaling limit of p(x,t) and 
their moments. We also obtain exact critical exponents for their moments at 
the critical point and show universality classes with respect to these critical 
exponents. 

1 Introduction 

In recent years quantum walks (QWs), as the quantum analog of the classical random 
walks (CRWs), have attracted great attention from mathematicians, computer scien- 
tists, physicists and engineers. Two forms of QWs, continuous-time QWs (CTQW) 
[7] and discrete-time QW (DTQW) [H EH EH S] , are widely studied. In this work, 
we restrict our discussion to DTQW. In this case, an extra "coin" degree of freedom 
is introduced into the system. Unlike the classical random walk, where the direction 
of the particle moves is determined by the outcome of a "coin flip", for quantum 
random walks both the "flip" of the coin and the conditional motion of the particle 
given by unitary transformations are needed. 

In 2003, Brun, Carteret and Ambainis jl] discussed the decoherent Hadamard 
walk on 1-dimensional integer lattice Z, and found the expressions for the first and 
second moments of the position and showed that in the long time limit the variance 
grows linearly with time with the diffusive character. However, they did not provide 
an exact expression for the higher order moments nor the limiting distribution, due to 
the complicated forms of the eigenvalues of the superoperator Ckk and the difficulty 
to evaluate the position probability analytically. In 2004, Grimmett, Janson and 
Scudo [8] obtained the scaling limit of quantum random walks without decoherence. 
Their scaling factor is | . Recently in [2] , the scaling limit of a quantum random walk 
with a Markov controlled coin process converges either with scaling factor | or 
depending on eigenvalue conditions of the walk operator. 

In this paper, we consider the model with decoherence operators given by [3]. We 
overcome the difficulties in there to obtain the scaling limit of decoherent quantum 
random walks, by analyzing P(u,t), the characteristic function (Fourier transforma- 
tion) of the position probability distribution p(x,t). Here p(x,t), x G Z, denotes 
the probability that the particle is found at position x at time t. It turns out that 
the convergence of P(-^,t) depends on the spectrum of the superoperator C kk . We 
show that if 1 is an eigenvalue of the superoperator with multiplicity one and there is 
no other eigenvalue whose modulus equals to 1, then P(^=,t) converges to a convex 
combination of normal distributions. In terms of position space, the rescaled prob- 
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ability mass function p t (x,t) = p(\/ix,t), x G Z/y/t, converges in distribution to a 
continuous convex combination of normal distributions. We give an necessary and 
sufficient condition for a U(2) decoherent quantum walk that satisfies the eigenvalue 
conditions. For the 0(2) quantum walks such as the Hadamard walk, and walks un- 
der real or complex rotations are discussed. An explicit limiting distribution formula 
is provided for these walks. 

Our article is organized as follows. In Section 2, we present basic concepts of 
quantum random walks. In Section 3, we present our main result about the limit 
of p t (x,t). In Section 4, we give examples that illustrate our results. We obtain 
the scaling limit and the exact critical exponents for their moments at the critical 
point. We show that the decoherent quantum random walks, with coin space unitary 
transformation U G 0(2), 9 ^ n — 0, 1, 2, 3, belong to the same universality class 
with respected to the critical exponents of all moments at their critical points. 

2 The unitary walk on the line and decoherence 

Consider a general quantum random walk on the 1-dimensional integer lattices Z. To 
be consistent, we adapt analogous definitions and notations as those outlined in [I]. 
We denote the state space by a Hilbert space H p ®H 2) where H p denotes the position 
space and H 2 denotes the coin space. The basis of the position space are \x >, where 
i62 and, the basis of the coin space are \R > and \L >. We will assume that the 
walk starts at the origin. The shift operators in H p are defined as follows 



where S~ and S + are unitary shift operators on the particle position. Let Pr, Pi be 
two orthogonal projections on the coin space H 2 spanned by \R > or \L >, where 
Pr + Pi = I. Let U be a unitary transformation on H 2 , that acts as the "flipping" 
of the coin. Then the evolution operator of the quantum random walk is given by 



S \x >= \x + 1 >, 




S \x >= \x — 1 >, 



(2.2) 



E = (S + ®P R + S~ ® P L )(I ® U) 



(2.3) 



The eigenvectors \k > of S , S + are 




(2.4) 



X 
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with eigenvalues 

S+\k >= e~ ik \k >, 

S~\k >=e ik \k > . (2.5) 
Therefore, in \k > basis, the evolution operator is 

E(\k > ®|$ >) = \k> ®{e~ ik P R + e ik P L )U\<$> >= \k > <g>C4|$ >, (2.6) 

where U k = (e~ tk Pn + e lk Pi)U is also a unitary operator. 

The decoherence on the coin space is defined as follows. Suppose before each 
unitary transformation acting on the coin, a measurement is performed on the coin. 
This measurement is given by a set of operators {A n } on H 2 which satisfy 

Y J KA n = L (2.7) 

n 

Through out this paper, we also assume that the measurement is unital, i.e., it 
satisfies 

^A n A* n = L (2.8) 

n 

After the measurement, a density operator \ on H 2 is transformed by 

X^X' = J2 A "* A n- ( 2 - 9 ) 

n 

The general density operator of quantum random walk is then given by 



P 



/dk f dk! , ,. . 

27 / 2^' >< 1 ® Xfcfc '' ( ^ 



where Xkk 1 £ L(H 2 ), and L(H 2 ) is a vector space of linear operators on H 2 . Then 
after one step of the evolution under coin space decoherence, the density operator 
becomes 

/dk C dk' 
\kxk'\®Y. U k A nXkk >A* n U* k , (2.11) 

Suppose the quantum walk starts at the state |0 > (g>|$o >i then the initial state 
is given by the density operator 



/dk f dk' . , , ,. , . . 

— —\kxk'\®\%><%\. (2.12) 



After t steps, the state evolves to 
dk [ dk' 



Pt 



j^j < ^\ k><k '\® E U k A nt ...U k A ni \%><%\A* ni U* k ,---A* n U* kl . 



m,...,nt 

(2.13) 

If we define the superoperator C kk > to be an operator which maps L^H?) to L(H2): 

L kk ,B = J2u k A n BA* n U* kf , V5 e L(H 2 ), (2.14) 

n 

then 

/dk f dk' ., . . i i . , 

27 / 2^ >< fc 1 ® £fefc ' l$0 >K $0 '- (2 ' 15) 
The probability of reaching a point x at time t is 

p(x,t) = Tr{{\x >< x\ ® I)pt} 

= 1 y d& y dfc' < x| A: >< fc'|a; > Tr{£^,|$ >< $ |} 

= j^jdkj dk'e^ k -^Tr{Cl k ,\<S> >< *o|}- (2-16) 

3 The limiting distributions of quantum walks with 
decoherence coin 

Let 

P(i/,i) =< e wx > t = J2e iux p(x,t) (3.1) 

ji- 
be the characteristic function of p(x,t). By the property of the 5 function 



2tt 
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and (12.161) . we have 

< e ivx >t = eiUX P( x > l ) 

X 

= E etra / S / d ^ X(k - k,) Tr{C\ k ,\% >< $o|} 

X 

= SyJ fE e "* M ^4^o >< $o|} 

= / § / ir 27r ^' - k - ^ Tr ^'\^ ><: 

= i_y cZfcTr{4 )fe+ J$ >< $o|}- (3.3) 
For any initial state O G L(H 2 ), the generating function of < e lvx > t is given by 

oo 

G(z,v) = J2 zt<eiVX> t ( 3 - 4 ) 



t=0 



1 /■ 00 

t=o 



(3.5) 



= J_ / d kTr{- 6}, (3.6) 

27T y 7 - ZL k ,k+u 

where |^| < 1 and O G L(H 2 ). Note that the generating function is well defined since 
the spectrum of C k>k+l/ is less than or equal to 1 by the following lemma. 

Lemma 3.1. Suppose U G ^(C) and the set of operators {A n } is unital. Let A be 
an eigenvalue of C kjk+U , then |A| < 1. 

Proof: Define \\C k , k \\ = sn^deL(H 2 ) with ||6|| 2 = Tr(6*6), then we 

have 

llA.fcll < 1 (see ref.[S]). 



\\C k , k+u 6\\ 2 = Tr((J2UkA n dAlUl +w y(J2u k A n dA* n U* k+v )) 

n n 

= T r (u k+u (Y,A n dA* n ru;u k (Y,A n dA* n )u* k+u ) 

n n 

= Tr(^2 A n OA* n )*(^2 A nOA* n )) 

n n 
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This shows that \\Cl 



k,k- 



\C 



k,k\ 



< 1 and therefore |A| < 1. 



It follows from the above lemma that Y2u^o( z ^-k,k+uYO converges in \z\ < 1 and 
thus I — zCk,k+u does not have any pole inside the disk \z\ < 1. 

It is well known that any O £ L(H 2 ) can be written as a linear combination of 
Pauli matrices: 

O = r I + T\U\ + r 2 a 2 + r 3 a 3 , (3.7) 
where a\ 2 3 = a x y z are usual Pauli matrices. Hence O can be represented by a 



column vector 



O 



\TZ) 



{3.1 



By Lemma I3TT1 and 1 — zCk,k+u £ L(L(H 2 )), we have 

G(z,u) 



< e lvx > t 



1 

2tH 



yt+l 



dz, 



'|z|=r<l z 

for some < r < 1. Let A be the matrix associated with 1 
to the Pauli matrices, then 



(3.9) 

zCk,k+u, with respect 



1-zC 



-O = A~ 



k,k+u 



fro\ 

n 

r 2 

w 



det(A) 



/A u 
A 12 
A 13 

\A U 



A 2 i 
A 22 
A 23 
A 2i 



A 31 
A 32 
A 33 
A 3A 



A 41 \ 
A 42 
A i3 
A A4 J 



fro\ 

n 

r 2 
\r 3 / 



where Ay is the cofactor of A. 



Note that Tr(o"j) = for i = 1,2,3, and Tr(ao) 
in (13. 6p only the first row action h(z, v) = A n r + 
Therefore 

2h(z,v 



dk- 



detA 



= 2. So when taking the trace 
^21^1 + A 3 ir 2 + A^Ta remains. 

-. (3.10) 



Let L = be the matrix representation of £k,k+u in terms of Pauli matrices. 

Then we have the following lemma. 

Lemma 3.2. Suppose U £ U 2 (C) and {A n } is unital. Then £k t k+v has the following 
representation 

( cosu x x x\ 
xxx 
xxx 
, % sin v x x x 



V 



7 
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Moreover, if ' v = 0, then has the following representation 



/l 0\ 

x x x 
x x x 
\0 x x x j 

Proof: Let U G U 2 {C). Then \detU\ = 1. Let detf/ = e^. We consider the 
normalized operator 



(3.11) 



Then W G SU 2 {C). By (|2~T4l) . £ fefe / is the same for [/ and Therefore, without 
loss of generality, we may assume that U G SU^ifC,) with the following form 



U 



a —j3 

(3 a 



e ik a -e ik /3 
e ik (3 



e tk a 



where a,P G C and \a\ 2 + \/3\ 2 = 1. Then 

£4 = 

and 

where i — 0,1, 2, 3. On the other hand, 
Hence 

hi(v) + Z 4 i(i/) 
/n(z/) - hiiy) 
+ Z31 0) 

ih\{v) + /3l(^) 



(3.12) 



Z u (i/) + hi{v) -ihi{v) + l&{v) 



k+v 



e lv 
e - * 



(3.13) 



S 



which implies that = cos(u), hii^) = hi {v) — and U^v) = isin(z/). In 

particular, if v = 0, then the first column of C k ,k is (1, 0, 0, 0) T . 

Next we finish our proof by showing lu(0) = 0, for i — 2,3, 4. Suppose 



UuKa^AlUt 



then 



hi{0) + hi{0) =t*(0), 
Zii(0)-Z4i(0)=ri(0), 



and hence 



Z«(0) = ^(ri(0) + r|(0)) = IrrCA.fctTi-i) = ^Trfo-i) = 0, (3.14) 

for i = 2,3, 4. The third equality in the above holds because £ k ,k preserves the trace. 

Theorem 3.1. Suppose U G ^(C), the set of operators {A n } is unital, 1 is an 
eigenvalue of £ k ,k with multiplicity one, and |A| < 1 for any other eigenvalue A of 
C Kk . Then 

lim P(-^=,t) = — / e- l ^ {0)u2 dk, Mv G [0,2tt], 

where Zo( z/ ) rooi of det(l — zC k;k+v ) = stzc/i i/iai Zo(0) = 1. 

By the Cramer-Levy Theorem (see e.g. Theorem 6.3.2 [6]), the above theorem 

implies that under the conditions, the distribution of the scaling limit is a continuous 

convex combination of normal distributions with variance Zq(0). In other words, if 

_z_ 



we define the rescaled probability mass function on -4= by 



p t (x,t) =p(Vix,t),x G — = , (3.15) 

then p t converges in distribution to the continuous convex combination of normal 
distributions whose density function is given by 

1 f 7 * 1 



2vr 



F(x) = — / e 2z o'w d/c, x G i?. 



o v^(0) 



9 



Proof: By ( 13. lip . U can be normalized to an SU2 operator W. Note that by 
(12.141) and (I2.16p . £kk+v and p(x,t) are the same for U and W. Therefore, without 
loss of generality, we may assume that U G ST^C). By (13 . 1 0[) and Cauchy's integral 
formula, 

2mJ lzl=r<1 z^ 1 2n J 2m J lzl=r<l z^detA 

Let zqM, ziM, Z2U'), Z3M be four roots of det A = 0. Then 



20 {") ' ziM ' 2 2 (f) ' z 3 (iy) 

are four eigenvalues of Ck,k+v By the assumptions that 1 is an eigenvalue of JLk,k 
with multiplicity one, and |A| < 1 for any other eigenvalue A of £k,k, we may make 
the following ordering 1 = |*„(0)| < |*i(0)| < |z 2 (0)| < |* 3 (0)|. Let I (z, v) = ^t^. 
By continuity of Zi(u),i = 0, 1, 2, 3, in both variables v and k, there exist a constant 
R > 1, and a small neighborhood V of v = such that |^o(i / )| < R < \zi(u)\ for any 
u E V and G [0, 27r], for all i = 1, 2, 3. By Cauchy's Residue Theorem, 

Note that there exists t such that E V for all 1/ G (0, 27r) if £ > to- Let ^(z, z/) = 
det A By compactness of {(z, v, k); \z\ = R,v G V,k G [0,27r]}, there exists a 
constant C such that | i- 'JV | < Con {(z, z/); |z| = R,v £ [0,27r],fc G [0,27r]}, for all 
t > tn. Therefore, 



Hence 



1 / 1 2li^ 1 
lim <b —-T-, frdz < lim — — C = 0. 



lim Res(l, z = 0) = — lim Res(l, z = zq{—=)). 

t— >oo «/£ 



For any fixed u, we have 



Let t —> 00, we have 



^ /, 2/i(l, 0) , / v \ 
lim ReslLz = z {— = )) = -3 lim z — f= 



-t-i 



dz 
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We claim that 

2/i(l,0) 



1; z'(0) = 0. (3.16) 



§f(M) 

By the Dominated Convergence Theorem, we have 

lim P(-^=,t) = — 1^ e^ z °(°> 2 dk. (3.17) 

yjt 27T J Q 

We will finish our proof by proving the claim (I3.16p . Let M denote the 3 by 3 
submatrix of L = (kj): 

M(u) 

Then the matrix 

A \ v=0 = ' h ~ zM(0) J • 
The cofactor A 2 \ = A 3 i = An = and An = det(/ 3 — zM(0)). By Lemma 
are eigenvalues of £k,k for i — 0, 1, 2, 3. Therefore U(o) ' T^Jyj are eigenvalues of 
M(0). Hence 

det(/ 3 - zM(0)) = (1 - -^-)(1 - -^-)(1 - -^-) (3.18) 

-21 (0) 2 2 (0) Z 3 (0) 

'2 - zi(0))(z - z 2 (0))(z - 2 3 (0)). (3.19) 



// 2 2 


h3 


I24 


= k2 








I43 


I44 


(\-z 








zi(0)z 2 (0)z 3 (0) 
On the other hand, 

TT^ 1 ' 0) = ^ L (n^ 1 ~ ^ l(0))(1 " Z2(0))(1 " ^ (0)) - (3 - 20) 
dz 2i(0)^(0)z 3 (0) 

Hence 

since O is a density operator. Next we will show that Zq(0) — 0. Since g(zo(u), v) = 0, 
we have 

_ dg(z (v),z) _ dg(z,v) i , dg(z,v) i , f . 

U- - - ^ U=, oH + ^ I^m^oW- 

11 



When v — 0, z (0) = 1, then the above equation becomes 



0=%^U + %^U4(0). (3.21) 

Consider the matrix A at z — 1, and note that /2i(^) = ^31 (^) = 0, we have 

g(l, is) = (l- lii(v))A n (v) + hx(v)A 21 (v) - l 31 (u)A 31 (u) + U^A^u) 
= (1 — cos(^))^4n(i/) + i sin vA A i(v). 

By Lemma [3. 2[ the cofactor A 41 (0) = 0. It follows that 

= °- (3.22) 

By ( I3.20p . ag g°'^ |z= Zo 7^ 0, since 2j(0) 7^ 1, for all i = 1, 2, 3, by the assumptions 
of the theorem. By (1X2TD and fl3T22"j) we have ^(0) = 0. 

Q.E.D 



4 Applications 

The assumptions that 1 is the largest eigenvalue of Ckk with algebraic multiplicity 1 
and that there is no other eigenvalues whose modulus equals to 1, are crucial in de- 
termining the convergence of P(^,t). In this section, we consider the measurements 
given by 

^0 = V 1 ~ P 1 ' 
A 1 = y /p\R><R\, 
A 2 = y/p\L >< L\. 

By analyzing the spectrum of Ckk, we obtain a necessary and sufficient conditions in 
which the assumptions of Theorem 13.11 are satisfied. Specific examples such as the 
Hadamard walk, walks under real and complex rotations are illustrated. For certain 
class of convergent quantum walks, explicit formulas are obtained for the limits of 
the characteristic functions of properly scaled p(x,t). In addition, we will also give a 
complete description of the behavior of those walks that do not satisfy the conditions. 
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Lemma 4.1. Let be a superoperator on the Hilbert space L(C 2 ), defined by 

2 

£ k , k <(d) = J2u k A n dA* n u* k „ 

n=0 

where U k and Uy are 2x2 unitary matrices and O G L(C 2 ). Then 

< jC k , k ,d,Ck,k>d ><< 6,6 > . 

In particular, < C k k >0, C^yO >=< 0,0 > if and only if the decoherence rate p = 
or 012 = o 2 i = 0. 

Part of the above lemma has been obtained by Liu and Petulante in [9], but 
our lemma extends the equality part of their lemma to a wider scope with more 
applications. 

Proof: The inequality follows from Lemma 1 in [9j. By the proof of Lemma 1 
in [9J, the equality holds if and only if (2p — p 2 )(o\ 2 + ^21) = 0- That is, p = or 
612 = 621 = 0. 

If p 7^ and 612 = 6 2 i = 0, then 1 may not be the unique eigenvalue of C k ,k with 
largest modulus. In this case, we have the following theorem. Let dim(X) denote the 
dimensions of the eigenspace associated with the eigenvalue A. 

Theorem 4.1. Let < p < 1. Let U G U(2). Suppose X is an eigenvalue of C k k, 
then we have 

( a ) E|A|=l^ m ( A ) < 2. 

(b) 1 is an eigenvalue of and dim(l) > 1. 

(c) If\X\ = 1, then X = 1 or — 1. 

(d) dim(l) = 2 if and only if u 12 = u 2 \ = and \uu\ = \u 2 2\ = 1. In this case its 
multiplicity is 2. 

(e) There exists eigenvalue X = — 1 if and only if u\\ = U22 = and \uu\ = |«2i| = 1- 
In this case its multiplicity is 1. 

Proof: By considering the normalized operator W as in (13. lip , and noting that 
the statements (a)-(e) do not depend on whether it is U or W, we mat assume without 
loss of generality that U is in SU 2 with the form 
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where a, (3 G C and \a\ 2 + \(3\ 2 = I. 

a) Let A be an eigenvalue of C kk with eigenvector 



021 °22 



Oil O12 



) 



.By Lemma 



13. 2\ Ck,k has the form 




(4.23) 



where 



2cos(2£;) J Re( / ga) + 2 sin(2£;) Im(pa), 
2cos{2k)Im{(3a) - 2 sin(2k) Re{f3a), 



(4.24) 
(4.25) 
(4.26) 
(4.27) 
(4.28) 



Z42 = 2qRe(a(3), 
Z 43 = 2iqlm(aj3), 
l M = \ a \ 2 - \(3\ 2 . 



Since |A| = 1, by Lemma we have 621 = 612 = 0. This implies that the dimension 
of the space spanned by the eigenspace for all eigenvalues with modulus 1 is at most 
2. Moreover, the intersections of eigenspace corresponding to different eigenvalues is 
{0}. Therefore a) holds. 

b) By (I4.23p . 1 is an eigenvalue of £k,k- Furthermore, (1,0,0,0) T is one of its 
eigenvectors. Therefore dim(l) > 1. 

c) Note that 

1 % 1 1 

O = -(du + 622)0-0 + -(612 - 621)0-1 + -(612 + 621)02 + -(6n - 622)03, (4.29) 

so if Lk t kO = AO, with |A| = 1, then by Lemma [4.11 




(4.33) 



(4.30) 



(4.31) 



(4.32) 
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(14.311) and (I4.32p can be written as the following matrix form 
'cos 2k — sin 2/c N 



sin 2k cos 2k 



(Re(3a, Im(3a) (6n — 622) = 0. 



(4.34) 



Note that the matrix on the left hand side of (I4.34p has determinant 1, for all k, 
hence invertible. Therefore, if 6 n — o 2 2 7^ 0, then \(3a\ = 0, or equivalently, a(3 = 0. 
Consequently, if /3 = 0, then A = 1. If a = 0, then A = — 1. On the other hand, if 
on — <5 2 2 = 0' then On + 6 2 2 7^ (otherwise O = 0). In this case, A = 1 by (14.30 p . 
So A = 1 or — 1. 

d) We first assume dim(l) = 2. This implies that there exists an eigenvector 
of the form (a, 0, 0, b) T , with b 7^ 0. Since (1, 0, 0, 0) T is an eigenvector, (0, 0, 0, 1) T 
is also an eigenvector. By (I4.29p . (1,0,0,0) T and (0, 0,0, 1) T are two eigenvectors 
corresponding to 6\\ — o 2 2 = and 6n + 622 = 0, respectively. Hence 
by (I4.33p . Therefore (3 = and |a| = 1. 

Conversely, if (3 = and \a\ = 1, then Ckk has the form 



1 



/I 








' 12 





U3 



0\ 




V 



(4.35) 



It follows from (14T261) and (14271) that / 42 = I43 = 0. Therefore, dim(l) > 2. So 
dim(l) = 2 by part a). Therefore, by (I4.35P again, the multiplicity of 1 is also two, 
otherwise dim(l) > 2 which contradicts a). 

e) If —1 is an eigenvalue, then its dimension must be 1 by part a). Suppose 
(a, 0, 0, 1) T is the associated eigenvector, then a = by (I4.30p . Therefore \a\ 2 — \(3\ 2 = 
-1 by (jQ3|) . This implies a = and \(3\ = 1. 

Conversely, if a = and \(3\ = 1, then & has the form 



/l 













123 
133 







-v 



(4.36) 



So dim(-l) = l by part a) and b). Therefore, by (14.361) again, the multiplicity of — 1 
is also one, otherwise dim(— 1) > 2 which contradicts a). 



Corollary 4.1. If U E 0(2), i.e. 
some 9 G [0,27r], then 



U 



cos 9 
sin# 



- sin 9 

cos 9 



or 



cos 9 
sin 9 



sin 9 
- cos 9 



for 
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a) 1 is an eigenvalue of Ck t k with multiplicity one, and |A| < 1 for any other 
eigenvalue holds if and only if 9 ^ ^ where n = 0, 1, 2, 3. 

b) If 9 = 0, it, 1 is an eigenvalue of Ck,k with multiplicity 2. 

c) If 9 = |, vCfc^ eigenvalues 1 and -1, each has multiplicity 1. 

Proof: Note that 9 = 0,ir if and only if w 12 = M21 = and |wn| = |?x 2 2 1 — 1) 
and 9 = |, ^ if and only if U12 = M21 = and |ttn| = \u22\ = 1- Therefore corollary 
follows. 

We are now ready to discuss examples according to different values of 9. 
In the Hadamard walk, the evolution operator is given by 

1 I ' p~ ik p~ ik \ 

(QQg — sin ~~ \ 
.4 4 tt )• By Theorem 14. 1[ the associated superoperator £k,k+u 

sin ^ — cos 2 J 

satisfies the assumptions in Theorem l3.lt and we have 

1 + q 2 + 2gcos2A; 



4(0) = 0; 4(0) 



where q = 1 — p. By Theorem I3.1[ we have 



q 2 



V 1 f 2lT 1 l+g 2 +2qcos2fc .,2 

lim P(-p,t) = — / e 5 l=? (ifc (4.38) 

A/t 27T J 

1 f 27T 1 l+g 2 +2 g cosfc ,,2 

= — / e 5 i=? dA;, (4.39) 
2tt Jo 

by change of variables and periodicity of cosine. 

In general, if U e 0(2) with det U = -1 and 9 ^ *f where n = 0, 1, 2, 3. Then 

y e Jfe sm(6 l ) -e tfe cos(6 l ) y ' v 7 

and 

/cos(z/) (1 — p)z sin(^) sin^) i sin(^) cos(26') \ 



C 



-(1 - p) cos(2£; + z/) cos(2#) (1 - p) sin(2A; + v) cos(2A; + 1/) sin(2#) 
-(1 -p)sin(2A; + cos (20) -(1 -p) cos(2A; + 1/) sin(2fc + v] |sm(20) 
\isin(i/) (1 — p) cos(z/) sin(2#) cos(i/) cos(20) J 
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which satisfies all the assumptions of Theorem 13.11 Direct computation gives 

, , , ii/ \ 1 + 2o cos(2A;) + q 2 0/ x 

z' (0 = 0; 4(0) = \ V cot g • 

1 — g z 



where q = 1 — p. Hence by Theorem I3.1[ 



- 1/ 1 /" l l+2qcos(2fc)+q 2 2( a\,,2 

lim P(-7=,t) = — / e~ 5 i=? cot dfc (4.41) 

1 /-2T _ii +2gcoB (fc )+g 2 a 2 



2vr 



e 



5 J=? cot (fl) " djfe, (4.42) 



o 



by change of variables and periodicity of cosine. 

Similar calculations also show that (14.421) holds for U £ 0(2) with det U = 1 and 
6 ^ ^,n = 0,l,2,3. 

The n-th moments M n of the limiting distribution can be calculated from C I4.42[) 
by using moment generation functions. Let 

2tt J 

Then 

oo 1 

V -i»M n i/ n = ^(i/) (4.43) 
^— ' n! 

n=0 

= E ^[(-^) n cot 2 "(g)^" — i -L Al + 2gcos(fc) + q 2 ) n dk (4.44) 

It follows that for U £ 0(2) with ^ , n = 0, 1, 2, 3, we have 

M 2n+1 = 0,n = 0,l,2,3,..., (4.45) 

and for even moments, 

M 2n = i-± cot 2 " W 2n(i _ g2)n - ^ (1 + 2gcos(2fc) + g 2 )^fc (4.46) 

^ cot 2 -(g) 2n(l _ g2)n - j( (e i2fc + g)"(e- i2fe + g)»dfc (4.47) 



^ cot2 " w nih^ S (') 2,2< """' " = °' 2 "" (4 ' 49) 
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«* Wtw; ^- / > > , e lU \ n ~ l , e" i( "g""' dfc (4.48 



Therefore we have 

(2n)l cot 2 9 
n\2 n 

where T n (q) is a polynomial of q of order 2n given by 



M 2n = ^)"T„(g), n = 0, 1, 2, ... (4.50) 



r -(9) = E(?)V- (4.5i; 



In particular, for Hadamard walk, the second moment of the limiting distribution 
is given in terms of 

T 1 (q) = l + q 2 . (4.52) 

This result agrees with the results given in jl]. 

Comparing to the well known 2n-th moment N2 n for the normal distribution with 
mean and variance a 2 = t i ot I , 

^ = ^(^)»,n= 0,1,2,..., (4.53) 

we see that the scaling limits of the decoherent quantum random walks are not 
normally distributed if q ^ 0. The deviation from the normal distribution gets larger 
as the even moments gets larger. However the deviations of the 2n-th moment are 
by the same factor T n (q) for all 9 ^ ^, where j = 0, 1, 2, 3. 

From (I4.50p . we also obtain the exact critical exponents for M 2n at p = 0: 

7 2n = hm = n. (4.54) 

p^o In p 

This result shows universality in which the critical exponents do not depend on 9 
as long as it converges. In other words, the coin-space decoherent quantum random 
walks, with coin space unitary transformation U G 0(2), 9 ^ n — 0,1, 2, 3, belong 
to the same universality class with respected to the critical exponents of all moments 
as p — > 0. 

Now we discuss the limiting behavior of the special cases. 

a) 9 = 0,7r. If the initial state is |0 > ®\R >, then E(\0 > ®\R >) = |1 > ®\R >. 
If the initial state is |0 > (g)|L >, then E(\0 > ®\R >) = -| - 1 > ®|L >. That is, 
the walk goes either left or right forever. Hence 



P(u,t) = e 



18 



depending on the initial conditions. In general, if the initial state is 0o = cr\R > 
+cl\L >, with \cr\ 2 + \cl\ 2 = 1, then 

P(j,t) = \c R \ 2 e™ + \c L \ 2 e-™ . 

b) 9 = f , ^f. If the initial state is |0 > ®\R >, then E(\0 > ®\R >) = -|1 > 
<8>\L >. If the initial state is |0 > (g>|L >, then E(\0 > ®\R >) = — | — 1 > <g>\R >. 
That is, the walk switches back and forth between two positions, which is trivial and 

lim P&,t) = l, 

t— >oo t 

for any k > 0. 



5 Concluding Remarks 

In this paper we consider coin space decoherent quantum random walks with coin 
space unitary transformation U. We prove that under the eigenvalue conditions, the 
scaling limit of the probability distribution converges in distribution to a continuous 
convex combination of normal distributions. An necessary and sufficient condition 
is obtained for U to satisfy the eigenvalue conditions. For U in 0(2), an exact form 
of the limiting distribution is given and the moments of all orders are obtained. For 
this case, the critical exponents are obtained and we show that all U with rotation 
angles 8 ^ n = 0,1,2,3, belong to the same universality class with respected 
to the critical exponents of all moments as p — > 0. Our analysis is based on the 
characteristic functions of the position distribution and the analysis of eigenvalues, 
Theorem 14. H and its corollary, which plays an important role in the applications of our 
main convergence theorem, Theorem 14. 11 We believe that a wider class of universality 
should hold for general quantum random walks in general d-dimensional lattices with 
general rotation U G U(n). For the future research, it would be very interesting to 
explore and classify their universality classes with respect to their critical points. On 
the other hand, we have fixed a measurement in our applications, while our general 
convergence theorem does not depend on the special form as that in our applications. 
An interesting problem would be to understand how the general measurements affect 
the limiting distributions and, especially their universality classes. 
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